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Abstract 

Purpose  -  Development  of  military  vehicles  capable  of  surviving  shallow-buried  explosive  blast  is 
seldom  done  using  full-scale  prototype  testing  because  of  the  associated  prohibitively  high  cost,  the 
destructive  nature  of  testing,  and  the  requirements  for  large-scale  experimental-test  facilities  and  a 
large  crew  of  engineers  committed  to  the  task.  Instead,  tests  of  small-scale  models  are  generally 
employed  and  the  model-based  results  are  scaled  up  to  the  full-size  vehicle.  In  these  scale-up  efforts, 
various  dimensional  analyses  are  used  whose  establishment  and  validation  requires  major 
experimental  testing  efforts  and  different-scale  models.  The  paper  aims  to  discuss  these  issues. 
Design/methodology/approach  -  In  the  present  work,  a  critical  assessment  is  carried  out  of  some 
of  the  most  important  past  efforts  aimed  at  developing  the  basic  dimensional  analysis  formulation  for 
the  problem  of  impulse  loading  experienced  by  target  structures  (e.g.  vehicle  hull)  due  to  detonation  of 
explosive  charges  buried  to  different  depths  in  sand/soil  (of  different  consistency,  porosity,  and 
saturation  levels). 

Findings  -  It  was  found  that  the  analysis  can  be  substantially  simplified  if  only  the  physical 
parameters  associated  with  first-order  effects  are  retained  and  if  some  of  the  sand/soil  parameters  are 
replaced  with  their  counterparts  which  better  reflect  the  role  of  soil  (via  the  effects  of  soil  compaction 
in  the  region  surrounding  the  explosive  and  via  the  effects  of  sand-overburden  stretching  and 
acceleration  before  the  associated  sand  bubble  bursts  and  venting  of  the  gaseous  detonation  products 
takes  place).  Once  the  dimensional  analysis  is  reformulated,  a  variety  of  experimental  results 
pertaining  to  the  total  blast  impulse  under  different  soil  conditions,  charge  configurations,  charge 
deployment  strategies,  and  vehicle  ground  clearances  are  used  to  establish  the  underlying  functional 
relations. 


Originality/value  -  The  present  work  clearly  established  that  due  to  the  non-dimensional  nature  of 
the  quantities  formulated,  the  established  relations  can  be  utilized  across  different  length  scales, 
i.e.  although  they  are  obtained  using  mainly  the  small-scale  model  results,  they  can  be  applied  at  the 
full  vehicle  length  scale. 

Keywords  Blast  impulse  loading,  Dimensional  analysis,  Shallow-buried  explosive  detonation 
Paper  type  Research  paper 

1.  Introduction 

Recent  efforts  of  the  US  Army  have  been  aimed  at  becoming  more  mobile,  deployable, 
and  sustainable  while  maintaining  or  surpassing  the  current  levels  of  lethality  and 
survivability.  Current  battlefield  vehicles  have  reached  in  excess  of  70  tons  due  to 
ever-increasing  lethality  of  ballistic  threats  which  hinders  their  ability  to  be  readily 
transported  and  sustained.  Therefore,  a  number  of  research  and  development 
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programs  are  under  way  to  engineer  light-weight,  highly  mobile,  transportable  and 
lethal  battlefield  vehicles  with  a  target  weight  under  20  tons.  To  attain  these  goals, 
significant  advances  are  needed  in  the  areas  of  light-weight  materials  and  structures  as 
well  as  in  the  areas  of  measurements  and  prediction  of  impulse  loads  resulting  from 
various  threats  (e.g.  LEDs,  buried  charges,  etc.).  The  present  work  focuses  on  the 
problem  of  quantifying  and  predicting  impulse  loads  resulting  from  buried  charges 
under  different  soil  conditions,  charge  configurations,  charge  deployment  strategies, 
and  vehicle  ground  clearances. 

The  main  issue  with  the  detonation-impulse  loading  measurement  efforts  is  that 
they  should  ideally  involve  full-scale  prototypes  of  the  vehicles,  which  is  quite 
expensive  and  requires  large-scale  test  facilities  and  a  large  crew  of  engineers 
committed  to  the  task.  On  the  other  hand,  it  is  quite  advantageous  to  carry  out  all 
experimental  testing  and  validation  using  small-scale  models  and  to  then  scale-up  the 
model  results  to  the  full-scale  prototype  using  the  appropriate  dimensional  analysis 
(defined  in  the  following  paragraph). 

The  main  objective  of  the  present  work  is  to  critically  assess  and  reformulate  the 
scaling  analysis  of  Wenzel  and  Hennessey  (1972)  for  military  vehicle-hull  floor  plates 
subjected  to  impulse  and  high-pressure  loads  associated  with  detonation  of  charges 
shallow-buried  in  soil  under  the  vehicle  hull.  The  dimensional  analysis  of  Wenzel  and 
Hennessey  (1972)  is  based  on  the  so-called  Buckingham  77 theorem  (Buckingham,  1914), 
a  procedure  commonly  used  in  the  situations  in  which  full-scale  structures/prototypes 
are  designed  and  sized  based  on  the  results  obtained  using  small-scale  model  tests. 
Within  the  critical  assessment  stage,  each  physical  parameter  (originally  identified  by 
Wenzel  and  Hennessey  (1972)  as  being  relevant  to  the  problem  at  hand)  is  re-examined 
with  respect  to  its  ability  to  provide  a  first-order  effect  to  the  impulse  loading  on  a  target 
plate  (a  surrogate  for  a  vehicle  hull).  In  the  dimensional-analysis  reformulation  stage,  an 
attempt  was  made  to  establish  a  functional  relationship  between  the  total 
charge-detonation  impulse  and  various  parameters  related  to  different  soil  conditions, 
charge  configurations,  charge  deployment  strategies,  and  vehicle  ground  clearance. 

The  main  role  of  the  military- vehicle  hull  floor  is  to  protect  the  vehicle  occupants  in 
the  case  of  shallow-buried  charge  detonation  underneath  the  vehicle  belly.  Complex 
interrelationships  generally  exist  between  the  type  and  the  extent  of  vehicle-occupants 
injuries  and  the  extent  of  vehicle  damage  resulting  from  under-the- vehicle 
shallow-buried  explosive  (SBE)  blast.  The  most  serious  vehicle-occupants  injuries  and 
their  fatalities  are  usually  incurred  as  a  consequence  of  hull-floor  rupture.  In  addition  to 
generating  floor-material  fragmentation,  floor  rupture  enables  the  ingression  of  flame 
and  toxic  (detonation-product)  gases  and  can  lead  to  on-board  fires  and  explosions. 
While  the  vehicle  under  shallow-buried  charge  attack  is  always  subjected  to  high  “G” 
forces,  this  aspect  of  the  threat  is  usually  considered  as  secondary  when  floor  rupture 
takes  place.  On  the  other  hand,  when  SBE  detonation  results  in  only  hull-floor  bulging, 
the  vehicle  is  then  propelled  upward  and  off  the  ground  to  an  extent  which  scales  directly 
with  size  of  the  explosive  and  inversely  with  vehicle  weight.  The  accompanying 
dynamic  shock  can  cause  vehicle  occupants  and  on-board  instruments  and  weapons  to 
be  thrown  about  within  the  vehicle  interior.  In  addition,  large  hull-floor  deflections  and 
deformations  can  cause  a  sequence  of  failures  throughout  the  entire  vehicle  through 
components/sub-system  connections  and  interfaces  (e.g.  via  fuel  lines  running  along  the 
floor,  floor-bolted  seats,  ammunition  storage  racks,  power-train  lines,  etc.). 


Traditionally,  the  floor-rupture  problem  is  solved  through  the  use  of  thicker  floor-plates, 
stronger,  tougher,  and  usually  heavier  material  and  through  the  utilization  of  applique 
armor.  This  approach  is  generally  associated  with  unnecessarily  oversized/overweight 
vehicle-hull  floors  and  is  being  abandoned  nowadays  due  to  aforementioned  military’s 
requirements  for  lighter  vehicles. 

It  is  clear  based  on  the  discussion  presented  in  the  previous  paragraph  that  in  order 
to  assist  vehicle  designers  in  the  development  of  lightweight,  high-survivability 
vehicle  hulls,  it  is  critical  to  be  able  to  accurately  measure  and  correlate  the 
charge-detonation  impulse  with  various  parameters  pertaining  to  different  soil 
conditions,  charge  configurations,  charge  deployment  strategies,  and  vehicle  ground 
clearances.  These  measurements  and  correlations  should  be  applicable  at  the  full 
vehicle-hull  length  scale  as  well  as  at  the  length  scale  of  smaller  models.  The  latter 
requirement  clearly  suggests  that  the  results  of  the  impulse  measurements  and  their 
correlations  with  various  parameters  related  to  different  soil  conditions,  charge 
configurations,  charge  deployment  strategies,  and  vehicle  ground  clearance  should  be 
defined  using  non-dimensional  quantities  which  transcend  different  length  scales.  In 
other  words,  dimensional  analysis  must  be  employed.  The  employment  of  dimensional 
analysis  in  correlating  impulse  loads  resulting  from  detonation  of  shallow-buried 
charges  with  various  parameters  pertaining  to  different  soil  conditions,  charge 
configurations,  charge  deployment  strategies,  and  vehicle  ground  clearance  is  the  main 
objective  of  the  present  work.  While  employing  the  dimensional  analysis  mentioned 
above,  three  sets  of  experimental  data  were  utilized  covering  a  large  range  of  target 
plate  sizes.  Details  regarding  the  test  setups  and  procedures  used  to  collect  these  data 
sets  are  presented  in  Section  3. 

It  is  hoped  that  application  of  the  proposed  dimensional  analysis  and  the 
associated  relations  between  the  physical  parameters  governing  the  problem  at  hand 
can  help  vehicle  designers  reduce  the  weight  of  the  vehicles  while  ensuring  the 
required  confidence  level  for  vehicle  survivability  and  survivability  of  vehicle 
occupants. 

The  organization  of  the  paper  is  as  follows:  overviews  of  the  Buckingham  77 
theorem  (Buckingham,  1914)  and  the  Wenzel  and  Hennessey  (1972)  dimensional 
analysis  are  presented  in  Sections  2.1  and  2.2,  respectively,  while  the  proposed 
reformulation  of  the  latter  analysis  is  discussed  in  Section  2.3.  A  brief  description  of  the 
experimental  setups  and  procedures  used  to  generate  the  data  needed  to  establish  the 
basic  functional  relationship  between  the  blast  impulse  and  various  soil  conditions, 
charge  configurations,  charge  deployment  strategies,  and  vehicle  ground  clearances  is 
provided  in  Section  3.  The  main  results  obtained  in  the  present  work  are  presented  and 
discussed  in  Section  4,  while  the  key  conclusions  resulted  from  the  present  study  are 
summarized  in  Section  5. 

2.  Background,  critical  assessment  and  reformulation 

2.1  Buckingham  T1  theorem 

As  discussed  earlier,  the  Buckingham  77  theorem  (also  known  as  “dimensional 
analysis”)  (Buckingham,  1914)  is  employed  in  the  present  work  in  order  to  address  the 
problem  of  scale-up  of  the  model  results  to  the  full-scale  prototype.  In  the  remainder  of 
this  section,  a  brief  overview  is  provided  of  dimensional  analysis  and  of  its  utility  and 
limitations. 
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Dimensional  analysis  is  essentially  a  technique/procedure  for  describing  the 
behavior  of  a  system  in  terms  of  a  number  of  dimensionless  parameters  (generally 
referred  to  as  the  77  terms).  This  procedure  typically  involves  the  following  main  steps: 

(a)  Identification  of  all  physical  parameters  which  govern  the  behavior  of  the  system 
under  investigation.  This  is  the  most  critical  step  in  the  dimensional  analysis 
since  it  requires  a  good  understanding  of  the  underlying  physics  of  the  problem 
at  hand  (even  in  the  case  when  the  functional  forms  of  the  governing  equations 
are  not  known). 

(b)  Specification  of  the  units  for  the  physical  parameters  defined  in  (a)  in  terms  of  the 
fundamental  physical  dimensions  (typically  force,  length,  and  time  or  mass, 
length,  and  time). 

(c)  Identification  of  a  subset  of  parameters  among  the  physical  parameters  defined 
in  (a)  with  the  number  of  parameters  in  this  subset  being  equal  to  the  number  of 
fundamental  physical  dimensions  defined  in  (b).  The  parameters  in  this  subset 
will  be  used  for  converting  the  remaining  parameters  in  (a)  into  a  set  of 
dimensionless  parameters  (i.e.  the  77  terms).  In  the  case  where  force/length/time 
or  mass/length/time  are  the  fundamental  physical  dimensions,  there  are  three 
physical  parameters  in  this  subset.  It  should  be  noted  that,  while  in  a  given 
problem  the  number  of  the  resulting  77  terms  is  uniquely  defined  (equal  to  the 
number  of  physical  parameters  in  (a)  minus  the  number  of  fundamental 
dimensions  in  (b)),  the  functional  form  of  the  77  terms  changes  with  the  choice  of 
the  physical  parameters  defining  the  subset  defined  earlier.  Furthermore,  it 
should  be  recognized  that  while  it  is  generally  advantageous  to  define  the  77 
terms  in  such  a  way  that  they  have  a  clear  physical  meaning,  regardless  of  the 
choice  of  the  77  terms,  the  final  implications/findings  offered  by  dimensional 
analysis  is  independent  of  the  expressions  chosen  for  the  77  terms  (provided,  all 
the  governing  physical  parameters  are  correctly  identified  in  (a)). 

(d)  Identification  of  the  functional  relationships  between  each  of  the  77  terms  and  all 
other  77  terms.  These  relations  are  generally  determined  using  experimental  or 
computational  approaches  or  the  combination  of  the  two  approaches. 

The  power  of  dimensional  analysis  is  that  the  functional  relations  offered  are 
generalized,  i.e.  the  effect  of  geometrical,  kinematic,  ambient,  loading  and 
response  scales  is  eliminated.  Consequently,  the  dimensional  analysis  predicts  that, 
if  all  the  pertinent  physical  parameters  were  correctly  identified  in  (a)  and  if  all  the  77 
terms  were  kept  invariant  (i.e.  unchanged)  between  a  small-scale  model  and  its 
full-scale  prototype,  then  the  results  obtained  using  a  small-scale  model  could 
be  directly  used  to  predict  the  response  of  the  corresponding  full-scale  counterpart.  The 
condition  under  which  the  77  terms  are  kept  invariant  is  generally  referred  to  as 
“Similitude  Theory”. 

2.2  Brief  overview  of  Wenzel  and  Hennessey  (1972)  work 

As  mentioned  earlier,  one  of  the  main  objectives  of  the  present  work  is  to  critically 
assess  and,  if  required,  reformulate  the  dimensional  analysis  presented  by  Wenzel  and 
Hennessey  (1972).  In  the  remainder  of  this  section,  a  brief  overview  is  provided  of  the 
dimensional  analysis  of  Wenzel  and  Hennessey  (1972). 


The  work  of  Wenzel  and  Hennessey  (1972)  deals  with  scaling  floor-plates  for  military 
vehicles  expected  to  be  targeted  by  SBE  blast.  The  main  objective  of  the  work  was  to 
develop  and  validate  a  dimensional  analysis  for  the  vehicle  floor  plate  sizing  problem,  so 
that  the  results  of  small-scale  replica-models  can  be  used  to  determine  the  minimal 
required  floor-plate  thickness  which  will  ensure  that,  under  the  specific  shallow-buried 
charge-detonation  scenario,  the  vehicle  floor  will  not  undergo  rupture  or  extensive 
dynamic  deflection.  The  term  “replica-model”  will  be  defined  in  greater  detail  later  in  this 
section.  Until  then,  this  term  will  be  used  to  denote  a  small-scale  model  which  is  similar 
(geometrically,  constitutively,  and  with  respect  to  extent  of  properly  scaled  loads  and 
structural  responses)  to  the  corresponding  full-scale  prototype. 

The  starting  point  in  the  work  of  Wenzel  and  Hennessey  (1972)  was  to  take 
advantage  of  the  fact  that  floor-plate  deflections  do  not  significantly  modify  the  (SBE 
detonation  induced)  dynamic  loading.  Toward  that  end,  the  problem  at  hand  was 
divided  (decoupled)  into  two  separate  problems: 

(1)  a  floor-plate  structural  response  problem;  and 

(2)  the  problem  of  analyzing  pressures  and  impulses  resulting  from  detonation  of 
an  explosive  (shallow  buried  in  soil  underneath  the  vehicle-hull  bottom). 

2.2.1  Structural  response.  In  their  work,  Wenzel  and  Hennessey  (1972)  carried  out  a 
careful  and  detailed  engineering  analysis  of  the  structural  response  of  a  thick  metal 
plate  (a  model  for  the  vehicle-hull  bottom)  supported  on  wooden  blocks  (to  obtain  the 
required  level  of  ground  clearance),  and  clamped  along  its  edges  (to  simulate  the  effect 
of  the  surrounding/supported  vehicle-frame  structure).  A  schematic  of  a  typical  test 
structure  used  by  Wenzel  and  Hennessey  (1972)  is  shown  in  Figure  1. 

Physical  parameters.  The  analysis  of  Wenzel  and  Hennessey  (1972),  mentioned 
above,  yielded  11  key  physical  parameters  which  govern  the  response  of  the  plate 
subjected  to  SBE  detonation  loading.  These  parameters  are  listed  and  explained  and 
their  fundamental  dimensions  are  identified  in  Table  I.  In  addition,  they  are  grouped 
in  accordance  to  their  type  (e.g.  input  blast  loading,  geometrical,  material,  response). 
The  parameters  which  define  the  loading  due  to  SBE  blast  are  the  (maximum) 
pressure,  P,  and  specific  impulse,  I.  Only  the  case  of  rectangular  target  plates 
was  considered  and,  hence,  the  physical  parameters  defining  plate  geometry  include 
thickness,  h,  and  the  two  edge  lengths,  Li  and  L2.  The  constitutive  and  inertial 


Dimensional 

analysis 


371 


Figure  1. 

A  schematic  of  the 
prototype  SBE-blast 
test  set-up  used  in  the 
work  of  Wenzel  and 
Hennessey  (1972) 
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Symbol  Parameter 

Fundamental 

dimensions 

Comments 

Input  blast-loading  parameters 

P  Blast  pressure 

F/L2 

For  the  rigid-plate  case,  the 

372 

I  Blast  specific  impulse 

FT/L2 

pressure  can  be  directly  related  to 
the  specific  impulse  and,  hence,  is 
an  irrelevant  parameter 

This  quantity  was  not  directly 

Plate  geometrical  parameters 
h  Plate  thickness 

L 

measured  in  the  present  work  and, 
hence,  was  eliminated.  Please  see 
Table  II  for  more  details 

Not  important  in  the  rigid-plate 

LhL2  Rectangular  target-plate  edge 

L 

case 

Control  the  fraction  of  detonation 

lengths 

Plate  material  parameters 
a  Plate  material  yield  strength 

F/L2 

momentum  transferred  to  the  plate. 
Li  is  equal  to  L2  in  the  case  of  a 
square  plate 

Not  important  in  the  rigid-plate 

ay  Array  of  other  strength  parameters 

_ 

case 

Not  important  in  the  rigid-plate 

for  plate  material 

e  Plate  material  equivalent  strain/ 

_ 

case 

Not  important  in  the  rigid-plate 

ductility 

p  Plate  material  density 

ft2/l 4 

case 

Not  important  in  the  rigid-plate 

Plate  response  parameters 

S  Plate  (bending)  deflection 

L 

case 

Not  applicable  in  the  rigid-plate 

a  Plate  acceleration 

L/T2 

case 

Related  to  the  blast  loading 

v  Initial  velocity  after  complete 

L/T 

pressure  and,  hence,  is  also 
irrelevant 

These  two  quantities  were  initially 

Table  I. 

momentum  transfer 

considered  as  they  are  directly 

Critical  assessment  of  key 

H  Plate  maximum  height 

L 

related  to  the  impulse  measurement 

physical  parameters 
controlling  structural 
response  of  the  plates 
originally  proposed  by 
Wenzel  and  Hennessey 

t  Post-SBE  detonation  time 

T 

methods  employed  in  the  present 
work.  However,  since  they  scale 
with  the  blast  impulse,  they  are  not 
used  in  this  analysis 

The  need  for  this  quantity  appears 

(1972)  and  the 
introduction  of  additional 
parameters 

unclear,  difficult  to  justify,  and 
hence,  is  not  utilized  in  the  present 
analysis 

response  of  the  target-plate  material  is  characterized  by  its  yield  strength,  several 
additional  stress-type  parameters  (e.g.  young’s  modulus,  ultimate  tensile  strength,  etc.), 
(Ti,  ductility,  e,  and  mass  density,  p.  The  deformation/motion  response  of  the  target 
plate  is  represented  using  the  maximum  flexural  displacement,  8  and  the  target  plate 
acceleration,  a.  Finally,  the  post-detonation  time,  t,  was  added  as  the  last  target  plate 
structural  response  parameter. 


77  terms.  Next,  Wenzel  and  Hennessey  (1972)  used  a  subset  of  parameters  consisting  of 
the  plate  material  yield  strength,  a,  the  plate  thickness,  h,  and  the  plate  material 
density,  p,  to  non-dimensionalize  the  remaining  physical  parameters  listed  above, 
creating  eight  77  terms. 

Scaling.  Finally,  within  the  structural  response  portion  of  the  problem,  Wenzel  and 
Hennessey  (1972)  established  the  basic  scaling  law  which  enables  the  use  of  small-scale 
models  to  obtain  results  applicable  to  a  full-scale  prototype.  According  to  this  scaling  law,  if 
the  linear  geometry  of  the  target-plate  model  is  scaled  by  a  factor  A  with  respect  to  the 
full-scale  structure  while  the  overall  shape  of  the  plate  and  the  plate  material  are  kept 
constant,  the  model  plate  deflection  will  scale  with  A,  the  maximum  blast  pressure  will 
remain  constant,  and  the  specific  impulse  will  scale  with  1/A.  In  accordance  with  the 
previous  discussion  regarding  similitude  theory,  when  the  model  is  constructed  in  such  a 
way  that  it  complies  with  this  scaling  law  (i.e.  when  the  model  is  a  so-called  “replica-model 
then  the  model-based  test  results  can  be  directly  applied  to  the  full-scale  prototype  if  each  77 
term  remains  unchanged  between  the  two  length  scales.  In  other  words,  both  the  model  and 
the  prototype  designs  are  associated  with  a  same  point  in  the  multi-dimensional  77  space 
and  the  knowledge  of  a  functional  relationship  between  the  77  terms  is  not  needed. 

2.2.2  SBE  detonation  induced  impulses  and  pressures.  Wenzel  and  Hennessey  (1972) 
also  carried  out  a  detailed  engineering  analysis  of  detonation  of  explosives 
shallow-buried  in  soil.  In  order  to  quantify  the  resulting  time-dependant  impulses 
and  pressures  (input  loading  parameters  in  the  structural  response  analysis),  they 
included  the  effect  of  the  explosive  shape,  size,  and  type,  as  well  as  properties  of  soil 
and  the  surrounding  (ambient)  air. 

SBE  related  parameters.  The  analysis  of  Wenzel  and  Hennessey  (1972),  mentioned 
above,  yielded  13  key  physical  parameters  which  govern  the  loadings  produced  by  a 
SBE  detonation.  These  parameters  along  with  their  fundamental  dimensions  are  listed 
in  Table  II.  In  addition,  they  are  grouped  in  accordance  to  their  type  (e.g.  charge, 
charge/plate  positioning,  ambient  air,  soil,  blast  loading,  etc.).  The  SBE  was  considered 
to  be  of  a  circular-disk  shape  and  to  have  a  size  defined  by  thickness,  t,  and  radius,  r.  The 
type  of  explosive  is  then  expressed  by  the  total  energy  released  during  SBE  detonation, 
E.  The  location  of  the  explosive  with  respect  to  the  model/prototype  is  defined  in  terms  of 
the  depth-of-burial,  d  (a  distance  between  the  top  surface  of  the  SBE  and  the  soil/air 
interface)  and  standoff  distance,  R  (distance  between  the  top  surface  of  the  charge  and 
the  bottom  of  the  test  plate.  To  account  for  the  effects  of  ah  blast,  three  ambient-air 
parameters  are  considered:  ambient  pressure,  Pot  sound  speed  in  ah,  a0,  and  a  ratio  of  the 
constant-pressure  and  constant- volume  specific  heats  in  air,  yB.  To  account  for  the  effect 
of  soil  (into  which  the  explosive  is  buried)  Wenzel  and  Hennessey  (1972),  assumed  that 
detonation  of  a  SBE  can  be  related  to  the  problem  of  landmine  detonation  and  the 
resulting  crater  formation.  Consequently,  they  assumed  that  two  soil-specific 
parameters  should  be  included:  the  (initial)  soil  density,  p,  and  its  seismic  velocity,  c. 
Two  parameters  are  introduced  associated  with  the  SBE-detonation  specific  impulse,  7, 
and  the  associated  (peak)  pressure,  P.  Lastly,  the  effects  of  gravity  are  included  through 
consideration  of  the  gravitational  acceleration,  g,  as  one  of  the  physical  parameters. 

77  terms.  Next,  Wenzel  and  Hennessey  (1972)  used  a  subset  of  parameters  consisting 
of  the  standoff  distance,  R,  the  explosive  internal  energy,  E,  and  the  seismic  velocity,  c, 
to  non-dimensionalize  the  remaining  physical  parameters  listed  above,  creating  ten  77 
terms. 
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Fundamental 

Symbol 

Parameter 

dimensions 

Comments 

SBE  blast-loading  output  parameters 

P 

Pressure  due  to  SBE  blast 

F/L2 

For  the  rigid-plate  case,  the  pressure  can 
be  directly  related  to  the  specific  impulse 
and,  hence,  is  an  irrelevant  parameter 

Isp 

Specific  impulse  due  to  SBE 

FT/L2 

This  quantity  was  not  directly  measured 

blast 

in  the  present  work  and,  hence,  was 
eliminated 

I 

Total  blast  impulse 

FT 

A  new  quantity  added  to  comply  with 
impulse  measurement  methods  employed 
in  the  present  work 

Charge  parameters 

E 

SBE  blast  released  energy 

FL 

Retained  in  the  present  analysis 

t 

Circular-disc  shape  SBE 
thickness 

L 

Is  not  expected  to  have  a  first-order  effect 

r 

Circular-disc  shape  charge 

L 

Controls  the  volume  of  soil  overburden 

radius 

ejected 

Vdet 

Explosive  detonation 

L/T 

A  new  quantity  added  which 

velocity 

controls  the  time  of  sand-overburden 
bubble  burst 

Charge/plate  positioning  parameters 

d 

SBE  depth  of  burial 

L 

Retained  in  the  present  analysis 

R 

SBE/plate  standoff  distance 

L 

Retained  in  the  present  analysis 

Ambient  air  properties 

Po 

Ambient  atmospheric 

F/L2 

Are  not  expected  to  have  a  first-order 

pressure 

effect  since  momentum  transfer  via  air 

a0 

Speed  of  sound  in  air 

L/T 

shock  is  relatively  small 

Yo 

Ratio  of  specific  heats  for  air 

- 

Soil  properties 

ft2/l 4 

P 

Density  of  soil 

Is  not  expected  to  have  a  first-order  effect 
because  the  soil  fracture  strain  (controls 
the  onset  of  bubble  burst  and  venting)  and 
irreversible-compaction  bulk  modulus 
(controls  the  absorption  of  detonation 
energy  by  the  soil  surrounding  the  charge) 
play  a  more  important  role 

c 

Seismic  velocity 

L/T 

Controls  the  volume/depth  of  compacted 
soil  surrounding  the  detonated  explosive 

Efr 

Soil  fracture  strain 

- 

A  new  quantity  added  which  controls  the 

Table  II. 

Critical  assessment  of  key 
physical  parameters 
defining  loading  from  a 
SBE  originally  proposed 
by  Wenzel  and 
Hennessey  (1972)  and  the 
introduction  of  additional 
parameters 


Kic  Soil  irreversible  compaction 
bulk  modulus 


onset  of  sand-overburden  bubble  burst 
and  venting.  For  the  same  type  of  soil,  it 
may  be  related  to  the  soil  porosity  and/or 
degree  of  saturation; 

F/L  2  A  new  quantity  added  which  controls 
the  amount  of  energy  absorbed  by  the 
soil  directly  surrounding  the  charge 
during  irreversible  compaction.  For  the 
same  type  of  soil,  it  may  be  related  to 
the  soil  porosity  and/or  degree  of 
saturation 

(continued) 


Symbol  Parameter 


Fundamental 

dimensions  Comments 


Dimensional 

analysis 


a  Soil  porosity 

P  Degree  of  saturation 

Gravity 

g  Gravitational  acceleration 


L/T2  Due  to  small  travel  times  and  extremely 


New  quantities  added  which  control  the 
maximum  degree  of  compaction 


375 


high  detonation  products/soil  ejecta 
velocities,  gravity  is  not  expected  to  have 
a  first-order  effect 


Table  II. 


Scaling.  As  in  the  case  of  the  structural  response  portion  of  the  problem,  Wenzel  and 
Hennessey  (1972)  established  the  basic  scaling  law  for  the  SBE  impulse-loading  portion 
of  the  problem  which  enables  the  use  of  small-scale  models  to  obtain  results  applicable 
to  a  full-scale  prototype.  According  to  this  scaling  law,  if  the  linear  dimensions  of  the 
charge  and  its  position  relative  to  the  plate  are  scaled  by  A  (the  same  scaling  factor 
used  in  the  structural  response  portion  of  the  problem),  while  the  overall  shape  of  the 
charge  as  well  as  the  ambient  air  and  soil  properties  are  kept  constant  and  gravity  is 
scaled  by  1/A,  then  the  maximum  blast  pressure  will  remain  constant  and  the  specific 
impulse  will  scale  with  1/A  (as  required  by  the  structural  response  portion  of  the 
problem).  In  accordance  with  the  previous  discussion  regarding  the  use  of  replica 
models  in  the  structural  response  portion  of  the  problem,  if  each  77  term  is  held 
invariant  between  the  two  length  scales,  then  the  model-based  results  can  be  directly 
applied  to  a  full-scale  case  and  no  knowledge  of  a  functional  relationship  between  the 
77  terms  is  needed. 

2.2.3  Critical  assessment  and  reformulation  of  Wenzel  and  Hennessey’s  (1972) 
dimensional  analysis.  The  main  objective  of  the  work  presented  in  this  section  is  to: 

•  critically  assess  the  dimensional  analysis  of  Wenzel  and  Hennessey  (1972)  in 
order  to  identify  the  physical  parameters  which  are  expected  to  have  a  first-order 
effect  on  the  impulse  loading  resulting  from  SBE  detonation;  and 

•  to  reformulate  the  analysis,  if  necessary,  in  order  to  introduce  77  terms  whose 
functional  forms  have  a  physical  basis  and,  consequently,  will  facilitate  the 
establishment  of  the  fundamental  relationships  between  them. 

A  brief  summary  of  the  results  obtained  in  this  section  is  provided  in  Tables  I  and  II 
while  a  more  detailed  description  of  the  work  carried  out  is  provided  in  the  following 
two  subsections. 

It  should  be  recognized  that  within  the  experimental  test  procedures  utilized  within 
the  present  work  to  establish  the  basic  functional  relations  between  77  terms,  the  target 
plate  used  was  thick  enough  that  the  plate  response  to  the  blast  loads  was  dominated 
by  the  rigid-body  motion.  In  other  words,  no  significant  flexion  or  damage  was 
incurred  by  the  target  plate.  Consequently,  the  identification  of  the  physical  parameters 
with  first-order  effects  was  carried  out  under  the  rigid  target-plate  assumption. 

Critical  assessment.  It  should  be  noted  that  the  parameters  identified  below  as  being 
irrelevant  or  unimportant  have  not  been  considered  further  in  the  present  analysis. 
Conversely,  the  parameters  identified  as  having  a  first-order  effect  have  been  retained: 
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Blast  pressure.  For  the  rigid-plate  case,  this  quantity  can  be  directly  related  to  the 
specific  impulse  and,  hence,  is  an  irrelevant  parameter. 

Blast  specific  impulse.  This  quantity  was  not  directly  measured  in  the  present 
work  and,  hence,  was  eliminated  and  replaced  with  the  total  impulse  delivered  to 
the  target  plate. 

Plate  thickness.  Not  important  under  the  rigid-plate  assumption. 

Rectangular  target-plate  edge  lengths.  These  parameters  have  a  first-order  effect 
since  they  control  the  fraction  of  detonation  momentum  transferred  to  the  plate. 
In  the  case  of  a  square  plate  (used  in  the  present  work),  these  lengths  are  equal. 
Plate  material  yield  strength.  Not  important  in  the  rigid-plate  case. 

Array  of  other  strength  parameters  for  plate  material.  Not  important  in  the 
rigid-plate  case. 

Plate  material  equivalent  strain/ductility.  Not  important  in  the  rigid-plate  case. 
SBE  plate  material  density.  Not  important  in  the  rigid-plate  case. 

Plate  (bending)  deflection.  Not  applicable  in  the  rigid-plate  case. 

Plate  acceleration.  Related  to  the  blast  loading  pressure  and,  hence,  is  also 
irrelevant. 

Post-SBE  detonation  time.  The  need  for  this  quantity  appears  unclear,  difficult  to 
justify,  and  hence,  is  not  utilized  in  the  present  analysis. 

SBE  blast  released  energy.  Has  a  first-order  effect  since  it  scales  directly  with  the 
total  blast  impulse. 

Circular-disc  shape  SBE  thickness.  Is  not  expected  to  have  a  first-order  effect. 

Circular-disc  shape  charge  radius.  Has  a  first-order  effect  since  it  controls  the 
volume  of  soil  overburden  ejected. 

SBE  depth  of  burial.  Has  a  first-order  effect  since  it  controls  the  volume  of  soil 
overburden  ejected. 

SBE/plate  standoff  distance.  Has  a  first-order  effect  since  it  controls  the  fraction 
of  ejected  soil  which  impacts  and  transfers  impulse  to  the  target  plate. 

Ambient  atmospheric  pressure.  Is  not  expected  to  have  a  first-order  effect  since 
momentum  transfer  via  air  shock  is  relatively  small. 

Speed  of  sound  in  air.  Is  not  expected  to  have  a  first-order  effect  since  momentum 
transfer  via  air  shock  is  relatively  small. 

Ratio  of  specific  heats  for  air.  Is  not  expected  to  have  a  first-order  effect  since 
momentum  transfer  via  air  shock  is  relatively  small. 

Density  of  soil.  Is  not  expected  to  have  a  first-order  effect  because  the  soil 
fracture  strain  (controls  the  onset  of  bubble  burst  and  venting)  and 
irreversible-compaction  bulk  modulus  (controls  the  absorption  of  detonation 
energy  by  the  soil  surrounding  the  charge)  play  a  more  important  role. 

Seismic  velocity.  Has  a  first-order  effect  since  it  controls  the  volume/depth  of 
compacted  soil  surrounding  the  detonated  explosive. 


•  Gravitational  acceleration.  Due  to  small  charge-to-target  travel  times  and 
extremely  high  detonation  products/soil  ejecta  velocities,  gravity  is  not  expected 
to  have  a  first-order  effect. 

Reformulation.  Based  on  the  discussion  presented  in  the  previous  subsection,  the 
following  six  parameters  originally  introduced  by  Wenzel  and  Hennessey  (1972)  are 
retained  in  the  present  analysis:  square  target-plate  edge  length,  SBE  blast  released  energy, 
circular-disc  shape  charge  radius,  SBE  depth  of  burial,  SBE/plate  standoff  distance,  and 
seismic  velocity.  To  this  list,  the  following  six  new  physical  parameters  are  added: 

(1)  Total  blast  impulse  (I).  A  new  quantity  added  to  comply  with  impulse 
measurement  methods  employed  in  the  present  work. 

(2)  Soil  fracture  strain  (s/r).  A  new  quantity  added  which  controls  the  onset  of 
sand-overburden  bubble  burst  and  venting.  For  the  same  type  of  soil,  it  may  be 
related  to  the  soil  porosity  and/or  degree  of  saturation. 

(3)  Soil  irreversible  compaction  bulk  modulus  (Kic).  A  new  quantity  added  which 
controls  the  amount  of  energy  absorbed  by  the  soil  directly  surrounding  the 
charge  during  irreversible  compaction.  For  the  same  type  of  soil,  it  may  be 
related  to  the  soil  porosity  and/or  degree  of  saturation. 

(4)  Soil  porosity  (a).  A  new  quantity  added  which,  together  with  the  degree  of 
saturation,  controls  the  maximum  degree  of  compaction. 

(5)  Degree  of  saturation  (jS).  A  new  quantity  added  which,  together  with  the  soil 
porosity,  controls  the  maximum  degree  of  compaction. 

(6)  Explosive  detonation  velocity  (v^t)-  A  new  quantity  added  which  controls  the 
time  of  sand-overburden  bubble  burst. 

It  should  be  noted  that  the  following  two  quantities  were  also  initially  considered: 

(1)  Initial  velocity  after  complete  momentum  transfer  (v).  This  quantity  was  initially 
considered  because  it  is  directly  related  to  the  impulse  measurement  methods 
employed  in  the  present  work.  However,  since  v  was  found  to  scale  directly  with 
the  blast  impulse,  it  is  not  used  in  the  present  analysis  (to  avoid  redundancy). 

(2)  Plate  maximum  height  ( H ).  This  quantity  was  also  initially  considered  because 
it  is  directly  related  to  the  impulse  measurement  methods  employed  in  the 
present  work.  However,  since  H  also  scales  with  the  blast  impulse,  it  is  not  used 
in  this  analysis  (to  avoid  redundancy). 

n  Terms.  The  12  physical  parameters  identified  in  the  previous  subsections  are  next 
normalized  using  the  following  subset  of  three  parameters  (seismic  velocity,  standoff 
distance,  and  total  explosive  energy)  to  form  a  set  of  nine  77  terms.  The  mathematical 
formulation  of  these  terms  is  provided  in  Table  III  while  the  corresponding  scaling  law 
is  given  in  Table  IV. 

3.  Experimental  procedures 

In  this  section,  a  brief  description  is  provided  of  the  experimental  setups  and 
procedures  used  by  McAndrew  (2009)  and  Cheeseman  (2009)  to  generate  blast 
total-impulse  data  (as  a  function  of  different  soil  conditions,  charge  configurations, 
charge  deployment  strategies,  and  target-plate  standoff  distances).  These  data  are 
analyzed  in  the  present  work. 
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17  Term 

Physical  meaning 

Similarity  type 

IT  =  L/R 

Square  plate  edge  length 

Geometrical 

n2  =  r/R 

Circular  disk-shaped  charge  radius 

Geometrical 

n3  =  d/R 

Scaled  depth  of  burial 

Geometrical 

U4  =  E/KicR3 

Non-dimensional  blast  energy  density 

Material 

n5  =  viel/c 

Non-dimensional  charge  detonation  velocity 

Material 

378 

n6  =  ic/E 

Non-dimensional  impulse  due  to  SBE  blast 

Loading 

Il7  —  a 

Soil  porosity 

Material 

Table  III. 

17  terms  controlling  total 

ns  =  p 

Degree  of  saturation 

Material 

Tig  =  Sfr 

Soil  fracture  strain 

Material 

impulse  delivered  to  the  Note:  Please  note  that  the  six  terms  appearing  in  this  table  were  obtained  by  proper 
target  plate  non-dimensionalization  of  the  quantities  appearing  in  Tables  I  and  II  using  E,  R,  and  c 


Symbol 

Parameter 

Scale  factor 

Original  parameters  of  Wenzel  and  Hennessey  (1972) 

L 

Square  target-plate  edge  length 

A 

r 

Circular-disc  shape  charge  radius 

A 

R 

SBE/plate  standoff  distance 

A 

d 

SBE  depth  of  burial 

A 

c 

Seismic  velocity 

1.0 

E 

SBE  blast  released  energy 

A3 

New  parameters  introduced  in  the  present  work 

I 

Total  blast  impulse 

A3 

Efr 

Soil  fracture  strain 

1.0 

K,c 

Soil  irreversible-compaction  bulk  modulus 

1.0 

Table  IV. 

a, 

Soil  porosity 

1.0 

Replica-to-prototype 

a 

Soil  degree  of  saturation 

1.0 

scaling  for  impulse 

Vdet 

Explosive  detonation  velocity 

1.0 

loading  on  a  target 

Parameters  introduced  in  the  present  work,  but  not  used 

plate  resulting  from 

V 

Initial  velocity  after  complete  momentum  transfer 

1.0 

blast  of  a  SBE 

H 

Plate  maximum  height 

A 

3.1  Vertical  impulse  measurement  facility 

The  Vertical  impulse  measurement  facility  (VIMF)  is  a  structural  mechanical  device 
located  at  the  Army  Research  Lab  at  Aberdeen  Proving  Ground,  MD  (ARL-APG)  that 
enables  direct  experimental  determination  of  the  imparted  blast-loading  impulse  via 
measurements  of  the  vertical  displacement  of  a  known  fixed-mass  vertical  guide  rail 
that  is  capped  with  a  1.524  m  X  1.524  m  X  0.101  m  steel  target  plate  which  serves  as  a 
momentum  trap  to  capture  the  blast  loading  of  the  buried  charge.  The  design  and 
operation  of  the  VIMF  has  been  described  in  detail  by  Gniazdowski  (2004),  Skaggs  et  al. 
(2007)  and  Taylor  et  al.  (2005)  and  will  be  only  briefly  discussed  here.  A  schematic 
drawing  of  the  VIMF  including  the  vertical  guide  rail  and  the  orientation  of  the  charge  in 
sand  at  different  levels  of  water  saturation  is  shown  in  Figure  2.  A  photograph  of  the 
VIMF  is  also  given  as  an  inset  in  Figure  2. 

Before  a  VIMF  test  is  conducted,  the  pit  underneath  the  target  plate  is  filled  with 
sand  of  a  given  saturation  level.  Typically,  commercially  available  (Quikrete)  sand  is 
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Figure  2. 

A  schematic  of  the 
vertical  impulse 
measurement  fixture 


used  which  consists  of  94.4  percent  sand,  0.3  percent  gravel,  and  5.3  percent  silt/clay. 
The  dry  density  of  this  sand  is  ca.  1.49  g/cm3  and  is  typically  associated  with  a  porosity 
level  of  0.35.  To  place  the  explosive  into  the  sand,  a  circular  hole  is  dug  in  the  middle  of 
the  pit.  After  placing  the  explosive,  the  excavated  sand  is  replaced  atop  the  charge  and 
compacted  to  form  a  sand  overburden  with  a  level  of  porosity  and  saturation  level 
comparable  to  that  of  the  surrounding  sand.  Before  the  charge  is  detonated,  a  probe  is 
inserted  into  the  soil  in  order  to  determine  the  saturation  level. 

To  create  the  required  water-saturated  sand  condition,  a  cylindrical  pit  3.65  m  in 
diameter  and  1.32  m  deep  is  first  constructed  in  the  soil  within  the  VIMF  test  area. 
Next,  approximately  14.2  m3  of  commercially  available  (Quikrete)  sand  is  placed  in  the 
pit.  The  sand  typically  consists  of  94.4  percent  sand,  0.3  percent  gravel,  and  5.3  percent 
silt/clay.  The  maximum  dry-sand  density  is  1.49  g/cm3  while  the  maximum  wet-sand 
density  is  1.91  g/cm3. 


3.2  Sandbox  test  facility 

The  VIMF  target  plate  of  dimensions  given  in  the  previous  section  represents  a  \  to 
5-model  of  a  typical  military  vehicle  hull.  To  generate  impulse  results  for  smaller-scale 
models,  the  sandbox  test  facility  at  ARL  is  used.  Essentially,  the  sandbox  facility  is  a 
smaller  and  less  elaborate  detonation-impulse  measurement  facility  than  the  VIMF 
with  typical  target  plate  dimensions  of  0.750  m  X  0.750  m  X  0.0127  m. 


4.  Results  and  discussion 

In  this  section,  a  summary  and  a  discussion  are  first  presented  of  two  sets  of  results.  The 
first  set  of  results  pertains  to  the  blast-impulse  measurements  carried  out  using  the  VIMF 
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or  the  sandbox  test  facilities.  The  second  set  of  results  pertains  to  the  soil  property  testing 
results  obtained  in  our  prior  work  (Grujicic  et  al.,  2006,  2007,  2008a,  b,  c,  d,  2009,  2010) 
which  were  used  to  assess  the  effect  of  soil  conditions  (primarily  porosity  and  the  degree 
of  saturation)  on  the  key  soil  properties  (the  seismic  velocity,  the  irreversible-compaction 
bulk  modulus,  and  the  strain  to  fracture).  Based  on  these  results  and  a  physical-level 
analysis  of  the  key  phenomena/processes  associated  with  detonation  of  SBE  s  and  with 
interactions  between  the  detonation  products  and  the  surrounding  soil  as  well  as  with 
interactions  between  the  detonation  products  and  soil  ejecta  with  the  target,  a  functional 
relation  is  proposed  for  the  non-dimensional  blast  impulse.  Finally,  a  multiple-regression 
linear  analysis  is  employed  to  parameterize  this  relation. 

4.1  Representation  and  overview  of  the  experimental  results 

In  this  section,  a  brief  summary  and  overview  of  the  experimental  data  generated  in 
McAndrew  (2009),  Cheeseman  (2009)  using  the  VIMF  and  sandbox  test  facilities  are 
given.  The  overview  encompassed  two  sets  of  VIMF  experimental  data  (McAndrew, 
2009;  Cheeseman,  2009)  and  one  set  of  sandbox  experimental  data  (McAndrew,  2009). 
Due  to  concerns  regarding  the  potential  misuse  of  the  data,  quantitative  aspects  of 
these  data  are  not  shown  in  their  entirety.  Instead,  the  data  are  described  in  more 
qualitative  and  physically-based  terms. 

The  first  set  of  the  VIMF -measured  blast-impulse  data  was  obtained  over  a  range  of 
soil  moistures  between  zero  to  20  percent,  a  range  of  standoff  distances  from  0.4  to 
1.45  m  and  a  range  of  charge  masses  from  0.5  to  6.0  kg.  All  remaining  charge 
parameters  such  as  charge  type  (TNT),  shape  (circular  disk  with  a  3:1 
diameter-to-height  ratio),  and  depth  of  burial  (0.102  m)  were  kept  constant. 

The  second  set  of  the  VIMF -measured  blast-impulse  data  was  obtained  by  varying, 
over  two  levels,  the  following  four  parameters,  where  the  two  levels  are  specified 
within  parentheses: 

(1)  charge  standoff  distance  (0.3556  and  0.4064  m); 

(2)  charge  mass  (2  and  4.45  kg); 

(3)  charge  depth  of  burial  (0.0762  and  0.1016  m);  and 

(4)  target  plate  area  (1.486  and  3.251  m2). 

As  far  as  the  soil  moisture  is  concerned,  it  was  varied  over  three  (0, 8  and  16  percent)  levels. 

The  sandbox  blast-impulse  data  were  obtained  by  varying  the  charge  mass  over  two 
levels  (0.25  and  0.5  kg)  and  the  standoff  distance  over  five  levels  (0.193,  0.2007,  0.2032, 
0.2286  and  0.3048  m),  while  the  moisture  content  was  held  constant  at  a  nominal  level  of 
5  percent.  Multiple  tests  were  carried  out  at  the  same  charge-mass/standoff-distance 
conditions  to  reveal  the  associated  impulse  variance. 

Examination  of  the  three  sets  of  impulse  data  revealed  that: 

•  As  expected,  detonation  impulse  increases  with  increasing  charge  mass. 

•  Detonation  impulse  decreases  significantly  with  increasing  standoff  distance  as 
a  result  of  the  detonation  energy  being  distributed  over  an  expanding 
hemispherical  shock  front,  resulting  in  a  lower  energy  per  unit  area  of  the  shock 
front  and,  hence,  a  lower  impulse  delivered  to  the  target  plate. 

•  The  impulse  transferred  to  the  target  plate  increases  with  increasing  moisture 
content  (e.g.  at  a  constant  charge  mass,  depth  of  burial  and  standoff  distance, 


increasing  the  moisture  content  from  0  to  20  percent  more  than  doubles  the  blast 
impulse).  This  finding  is  related  to  the  fact  that  at  higher  moisture  contents,  the 
volume/mass  of  the  ejected  soil  is  larger,  the  soil  can  undergo  less  irreversible 
compaction  (and,  hence,  absorb  less  detonation  energy)  and  possesses  a  higher 
fracture  strain  (ductility),  allowing  the  sand  overburden  to  be  accelerated  to 
higher  velocities  by  the  expanding  gaseous  detonation  products  (i.e.  higher  soil 
ductility  delays  the  onset  of  the  overburden  burst). 

•  As  expected,  as  the  size  of  the  target  plate  increases,  the  impulse  captured  by  the 
target  plate  also  increases. 

4.2  Derivation  of  the  functional  relationships  between  II  terms 
In  this  section,  an  attempt  is  made  to  construct  and  parameterize  a  functional 
relationship  between  the  FI  terms  governing  the  problem  of  detonation  of  SBE  s  and  the 
interactions  between  the  detonation  products,  soil  ejecta  and  the  target  plate  (under  the 
rigid- target  condition).  While  it  is  common  practice  in  the  field  of  dimensional  analysis  to 
construct  these  relations  using  a  monomial  functional  form  applied  directly  to  the  77 
terms,  such  an  approach  typically  obscures  the  essential  physics  of  the  attendant 
phenomena/processes.  Instead,  physical  arguments  and  insights  will  be  used  in 
constructing  this  relationship  in  the  present  work. 

The  basic  relationship  to  be  derived  will  have  the  non-dimensional  impulse 
(represented  as  776  on  the  left-hand  side)  while  the  remaining  77  terms  listed  in  Table  III 
will  appear  on  the  right-hand  side.  Within  the  present  work,  the  terms  on  the  right-hand 
side  of  this  relation  will  be  organized  in  four  groups  with  each  group  accounting  for  a 
specific  physical  phenomenon.  The  four  physical  phenomena  identified  as  having  a 
first-order  effect  are: 

(1)  the  relative  area  of  the  target  plate  to  that  of  the  charge; 

(2)  the  volume  of  overburden  sand  ejected; 

(3)  the  energy  released  during  detonation  of  a  SBE;  and 

(4)  the  energy  absorbed  by  the  sand  surrounding  the  charge  during  irreversible 
compaction  taking  place  prior  to  overburden  bubble  burst  and  venting. 

The  relative  area  of  the  target  plate  to  that  of  the  charge.  This  quantity  controls  the 
fraction  of  soil  ejecta  which  impinges  upon  the  target  plate.  Considering  the  fact  that 
typically  in  excess  of  70  percent  of  the  total  impulse  transferred  to  the  plate  is  the  result 
of  soil  ejecta  impingement  upon  the  plate,  this  quantity  is  identified  as  having  a 
first-order  effect  on  the  total  impulse  transferred.  To  account  for  this  effect,  the  following 
simple  model  for  the  first  term  on  the  right-hand  side  of  the  776  relation  is  proposed:  1  — 
exp(-  5.8(77Mi  -  Ilf}).  The  term  IIi  -  772  represents  the  target-plate  overhang  (i.e.  the 
extent  by  which  the  target  half-edge  length  exceeds  the  charge  radius)  divided  by  the 
standoff  distance  and  the  coefficient  5.8  accounts  for  the  fact  that  when 
the  tan  ~  ’(77,  -  772)  is  about  50°,  the  target  plate  captures  nearly  all  (i.e.  99.9  percent) 
of  the  blast  impulse  carried  by  soil  ejecta  (Grujicic  et  al.,  2006). 

The  volume  of  overburden  sand  ejected.  As  discussed  earlier,  the  majority  of  the 
impulse  transferred  to  the  target  plate  takes  place  through  impingement  of  soil  ejecta 
upon  the  target-plate’s  strike/lower  surface.  Consequently,  the  total  impulse  transferred 
to  the  plate  should  scale  with  the  total  volume  of  soil  ejected  where  overburden  soil 
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constitutes  the  majority  of  this  volume.  For  the  circular-disc  shape  of  the  buried  charge, 
the  following  relation  is  proposed  for  the  second  term  on  the  right-hand  side  of  the  776 
relation:  773  III  exp(l  ,79/7g).  The  term  /73  77f  represents  the  volume  of  the  soil 
overburden  divided  by  the  third  power  of  the  standoff  distance,  the  term  exp(}.79TIs) 
accounts  for  the  effect  of  moisture  content  on  the  volume  of  ejected  soil  and  the 
coefficient  1.79  accounts  for  the  fact  that  the  volume  of  ejected  fully-saturated  (/3  =  1.0) 
soil  is  typically  about  six  times  larger  than  the  corresponding  ejected  volume  of  dry 
{fi  =  0.0)  soil  under  identical  explosive-charge  deploy ment/detonation  conditions 
(Grujicic  etaL,  2006). 

The  energy  released  during  detonation  of  a  SBE.  To  construct  a  term  pertaining  to 
the  energy  released  during  detonation  of  a  SBE,  the  approach  borrowed  from 
numerous  investigations  dealing  with  air  blast  (Pehrson  and  Bannister,  1997)  is 
utilized  in  the  present  work.  Within  this  approach,  the  total  energy  released  by 
detonation  is  divided  by  a  third  power  of  the  standoff  distance  in  order  to  obtain  a 
quantity  which  scales  with  the  associated  volumetric  energy  density  within 
the  hemispherical  shock  front  at  the  instant  when  the  expanding  outer  surface  of 
the  hemisphere  reaches  the  target  plate.  The  resulting  E/R  J  term  has  units  of  stress 
and,  hence,  has  to  be  non-dimensionalized  by  dividing  it  with  a  stress-based  quantity 
(the  irreversible  compaction  bulk  modulus  of  the  soil,  Klc,  in  the  present  work).  This 
non-dimensionalization  is  proposed  in  order  to,  at  least  partly,  include  the  effect  of 
energy  absorption  by  the  soil  surrounding  the  charge  during  irreversible  compaction. 
Thus,  the  third  term  on  the  right-hand  side  of  the  776  relation  can  be  written  as  774. 

The  energy  absorbed  by  the  surrounding  sand.  It  should  first  be  recognized  that  a 
portion  of  the  energy  absorbed  by  the  sand  surrounding  the  charge  during  irreversible 
compaction  taking  place  prior  to  overburden  bubble  burst  and  venting,  Eabsorbed,  was 
already  included  in  the  third  term  through  the  use  of  the  irreversible  compaction  bulk 
modulus.  In  addition  to  KIC,  the  other  terms  affecting  the  energy  absorbed  by  the  soil  are: 

•  the  square  of  the  soil  volumetric  strain  at  full  compaction,  [a(l  -  /3)]2;  and 

•  the  volume  of  the  fully-compacted  soil,  which  scales  with  the  cube  of  the  product 
of  c  and  tbursl  where  tburst  is  the  post-detonation  time  at  which  the 
sand-overburden  bubble  bursts  and  venting  takes  place. 

Using  a  simple  kinematic  analysis  of  the  sand  overburden-bubble  evolution  (driven  at 
a  velocity  Vdet  by  the  expanding  gaseous  detonation  products),  this  post-detonation 
time,  tburst,  was  found  to  be  equal  to  e°r5  •  r/vdet.  The  fourth  term  on  the  right-hand  side 
of  the  776  equation  is  then  defined  as: 

1  -  E absorbed /E  =  1  -  0.5 -Kk  ■  [a(l  -  /3)]2  ■  4/3 '77 ■  (c  ■  tf/E 

=  1  -  2/3  •  77 •  Kic  ■  Ml  -  p)]2  ■  c3  •  4/2  •  r3/[*4  •  E] 

=  1  -  2/3 ■  77-  n2  ■  (1  -  77g)2  •  if2  ■  nf/[if5  •  tt4] 

Proposed  expression  for  the  non-dimensional  impulse,  n6.  Based  on  the  aforementioned 
consideration  of  the  four  phenomena  identified  as  having  a  first-order  effect,  the 
following  general  expression  is  originally  considered  for  the  non-dimensionalized 
impulse  776: 


n6  =  ci(i  -  e_5-8(ni_n2))C2(n3n^L79Il8)C3(n4)C4 

X  (1  -  y  n?(i  -  n8)2n8nf /(n8n4)  (1) 

As  can  be  seen  in  this  relation,  the  monomial  form  is  still  used  but  it  is  applied  to  the 
four  physical  terms  identified  above  rather  than  to  the  individual  77  terms.  After 
carrying  out  a  preliminary  evaluation  of  the  five  coefficients,  appearing  in  equation  (1), 
Ci  -  c5,  it  was  found  that  the  number  of  independent  coefficients  can  be  reduced  to 
three,  i.e.  c4,  c2  ~  2-c4,  c3,  c4,  c5  ~  c4.  The  three  independent  coefficients  are  next 
renamed  as  follows:  a\  =  c4 ,  a2  =  c3,  a3  =  c4,  and  the  776  relation  re-written  as: 


rig  —  «i  n3rv 


(i 


2  1  79n8\ffl2 


~5.8(n1-n2)'\2 


)  (n4)  i 


-n  (i 


n8)2n8nf/(n2n4) 


(2) 


4.3  Functional  relationships  pertaining  to  the  effect  of  soil  type/condition 
Equation  (2)  contains  a  number  of  77 terms  which  are  directly  or  indirectly  affected  by  the 
type  and  physical  conditions  (e.g.  porosity,  degree  of  saturation,  average 
particle  size,  particle  size  distribution,  etc.)  of  the  soil.  In  the  present  work,  the  same 
soil  type  (Quikrete)  was  used  so  that  the  primary  parameters  of  the  soil  are  the  porosity 
and  the  degree  of  saturation.  Before  a  multiple-regression  analysis  can  be  applied  to 
equation  (2)  in  order  to  complete  evaluation  of  the  three  776-functional  coefficients 
( d\  -  fl3),  one  must  provide/determine  functional  relationships  for  the  soil-related 
parameters/77- terms  appearing  in  equation  (2).  This  was  done  here  by  utilizing  the 
results  of  our  CU-ARL  soil-model  development  work  reported  in  a  series  of  publications 
(Grujicic  et  al. ,  2006, 2007, 2008a,  b,  c,  d,  2009, 2010).  Based  on  the  results  obtained  in  that 
work  and  using  a  simple  strength-of-materials  type  of  analysis,  the  following  results  are 
obtained. 

Seismic  velocity.  The  effect  of  soil  porosity  (a)  on  the  speed  of  sound  in  soil  at  two 
levels  of  saturation  (/ 3  =  0.0,  dry  soil  and  (3  =  1.0,  saturated  soil)  is  shown  in  Figure  3. 
A  simple  least-squares  based  curve-fitting  procedure  was  applied  to  the  results 
obtained  by  Grujicic  et  al.  (2006,  2007,  2008a,  b,  c,  d,  2009,  2010)  which  yielded  the 
following  functional  relationship: 


C  Cfuii  .  compaction  ( 1  u) 


4.1 — 0.3/3 


(3) 


where  the  sound  speed  at  full  compaction  of  the  soil,  CfUu.COmpacUom  is  set  to  4,620  m/s. 

Irreversible-compaction  bulk  modulus.  Kic  is  defined  as  a  negative  rate  of  change  of 
pressure  with  a  change  in  the  volumetric  strain  when  the  volumetric  strain  is 
associated  with  irreversible  compaction/densification  of  the  soil.  KIC  is  assumed  in  the 
present  work  to  be  given  as  a  degree-of-saturation  (/3)  weighted  average  of  the 
corresponding  quantities  in  dry  soil,  Kic  dry  and  Kic  saturated.  The  latter  two  quantities 
are  found  to  depend  on  a  and  /fin  the  following  ways: 


Kic.dry  F-icdryfu]i—compacfioneXp^ 


(4) 
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Figure  3. 

The  effect  of  soil  porosity 
(a)  and  the  degree  of 
saturation  (j3)  on  the 
seismic  velocity 


where  K ic.dry, full-compactions  (  —  17.5  GPa),  Kgio2  ( —  57  GPa),  Kh2o  (  —  2.26  GPa),  are, 
respectively,  the  irreversible-compaction  bulk  modulus  for  fully  compacted  dry  sand, 
elastic  compression  bulk  modulus  of  Si02  and  elastic-compression  bulk  modulus  of 
water. 

The  soil  irreversible-compaction  bulk  modulus  at  an  intermediate  level  of  saturation 
is  then  defined  as: 


Kk  —  (1  (3)K  u:.dry  +  /Iff  tr. saturated  (6) 

The  effect  of  porosity  and  the  degree  of  saturation  on  the  irreversible-compaction  bulk 
modulus  as  defined  by  equation  (6)  is  shown  in  Figure  4. 

Fracture  strain.  Our  prior  work  (Grujicic  et  al,  2006,  2007,  2008a,  b,  c,  d,  2009,  2010) 
established  that  while  the  degree  of  soil  compaction  (i.e.  porosity)  affects  soil 
strength  (including  its  fracture  strength),  it  has  no  significant  effect  on  the  soil 
ductility  (as  measured  by  the  strain  to  fracture).  On  the  other  hand,  through 
capillary-enhanced  inter-particle  cohesion,  the  degree  of  saturation  has  a  major  effect  on 
the  soil  fracture  strain,  Figure  5.  Using  the  results  obtained  in  our  prior  work 
(Grujicic  et  al.,  2006,  2007,  2008a,  b,  c,  d,  2009,  2010),  fracture  strain,  s/r,  is  found  to 
depend  in  the  following  way  on  the  degree  of  saturation: 


Efr  Sfr  .saturated  saturated  &fr,dry)@ 


-4.5/3 


(7) 


where,  Sfr,dry  (  =  0-06)  and  efKsatumted  (  =  0.11),  are  the  dry  and  saturated  soil  fracture 
strains. 


4.4  Linear  regression  analysis 

To  determine  the  three  776-functional  coefficients,  a\  -  a3,  a  logarithmic  mathematical 
operation  is  applied  to  both  sides  of  equation  (2)  and  a  simple  multi-regression  linear 
analysis  is  employed  to  the  transformed  equation  (2)  in  which  both  sides  are  evaluated 
using  the  aforementioned  three  sets  of  experimental  data  (McAndrew,  2009; 


Note:  Please  note  that  while  soil  porosity  (a)  affects  soil  strength, 
it  has  no  first-order  effect  on  the  soil  ductility,  i.e.  the  fracture 
strain 


Cheeseman,  2009).  This  procedure  yielded  the  values  for  the  three  776  functional 
coefficients,  a\  -  a 3.  However,  due  to  previously  mentioned  concerns  regarding  the 
misuse  of  the  results  generated  in  this  work,  the  values  for  these  coefficients  are  not 
reported  here. 
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Figure  4. 

The  effect  of  soil  porosity 
(a)  and  the  degree  of 
saturation  (/3)  on  the 
irreversible-compaction 
bulk  modulus 


Figure  5. 

The  effect  of  the  degree  of 
saturation  (/3)  of  soil  on  the 
soil  fracture  strain 
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An  example  of  the  results  obtained  in  the  multi-regression  linear  analysis  is  shown 
in  Figure  6  in  which  a  plot  is  given  showing  the  correlation  between  measured  776 
values  and  their  computed  counterparts.  The  goodness-of-fit  of  the  measured  impulses 
by  the  proposed  functional  relation,  equation  (2),  is  quantified  by  the  Pearson 
correlation  coefficient,  Rcorr,  as  calculated  using  the  formula: 


Rn 


1  _ -A  / Vi 

w-i-AfV  s*  ;V  sy  ) 


(8) 


where  :r  =  lTgp  andy  =  If™1p  are  used  to  denote,  respectively,  the  measured  (McAndrew, 
2009;  Cheeseman,  2009)  and  the  computed  (equation  (2))  non-dimensional  impulse  values,  n 
represents  the  total  number  of  measured/computed  data  points,  an  over-bar  designates  the 
mean  value,  and  5  represents  the  corresponding  standard  deviation. 

The  computed  value  of  the  correlation  coefficient  ( Rcorr  =  0.90)  appears  to  be 
relatively  small,  suggesting  that  perhaps  the  proposed  impulse  function,  equation  (2),  is 
not  the  most  appropriate.  However,  before  reaching  this  conclusion,  one  must  account 
also  for  a  relatively  large  range  of  values  of  the  measured  blast  impulse  under  nearly 
identical  test  conditions  (i.e.  for  the  large  blast-impulse  variance).  To  help  clarify  this 
point,  a  few  horizontal  error  bars  are  drawn  in  Figure  6  to  denote  the  typical  range  of 
experimentally-measured  776  values  obtained  under  almost  identical  test  conditions. 
Since  in  most  cases,  the  rF£'np  vs  lT6xp  correlation  line  intersects  these  error  bars,  one 
should  conclude  that  despite  a  relatively  low  value  of  the  correlation  coefficient,  the 
proposed  blast-impulse  function,  equation  (2),  appears  appropriate.  This  conjecture  is 
further  supported  by  a  finding  that  if  the  individual  77gXp  values  associated  with  the 
nearly-identical  test  conditions  are  replaced  with  their  respective  mean  values,  the  value 
of  the  correlation  coefficient  increases  to  Rcorr  =  0.94.  It  should  also  be  pointed  out  that 


Figure  6. 

A  correlation 
analysis  between  the 
measured  non-dimensional 
blast  impulse  data  and 
their  corresponding 
computed  values  obtained 
using  the  dimensional 
analysis  developed  in  the 
present  work 


Measured  Non-Dimensional  Impulse 


some  of  the  VIMF  data  were  obtained  for  the  case  of  a  V-shaped  target  plate.  Since  the 
analysis  carried  out  in  the  present  work  assumed  that  the  target  plate  is  flat,  the 
experimental  results  in  question  are  not  expected  to  correlate  well  with  the  proposed  IJ6 
relation.  This  expectation  was  confirmed  by  the  finding  that  the  correlation  coefficient 
increased  from  0.90  to  0.92  when  the  experimental  data  corresponding  to  the  V-shaped 
target  plate  were  not  considered  in  the  multiple-regression  linear  analysis. 

To  quantitatively  validate  the  scaling  aspect  of  the  dimensional  analysis  carried  out 
in  the  present  work,  a  multiple-regression  linear  analysis  case  was  also  considered  in 
which  only  the  VIMF  data  were  used.  This  procedure  yielded  a  correlation  coefficient 
Rcorr  =  0.92.  Since  this  value  is  somewhat  higher  than  the  value  Rcorr  =  0.90  (based  on 
the  consideration  of  both  the  VIMF  and  sandbox  results),  one  may  question  the  validity 
of  the  present  scaling  analysis.  However,  one  must  also  take  into  consideration  the  fact 
that  in  most  cases,  the  “rigid-plate”  condition  was  not  satisfied  for  the  sandbox  facility 
experiments.  While  target  plate  deflections  are  typically  assumed  to  not  alter  the 
magnitude  of  the  blast  load  (W enzel  and  Hennessey,  1972),  they  may  seriously  affect 
the  accuracy  with  which  the  blast  impulse  is  measured.  In  fact,  it  was  found  that 
leaving  out  a  few  of  the  “outlier”  sandbox  data  points  can  make  the  difference  between 
the  two  values  of  the  correlation  coefficient  smaller  than  0.005. 
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4.5  Application  of  the  proposed  n6  relation 

In  the  previous  section,  it  was  established  that  the  proposed  776  relation  can  reasonably 
well  account  for  the  observed  effects  of  various  soil  conditions,  charge  configurations, 
charge  deployment  strategies  and  vehicle  ground  clearances  on  the  measured  blast 
impulse  values.  In  this  section,  a  brief  description  is  provided  regarding  the  application 
of  the  776  relation,  equation  (2).  In  assessing  impulse  loads  experienced  by  a  target 
structure  during  a  fully-described  landmine  blast  scenario,  the  following  procedure 
should  be  utilized: 

•  All  parameters  listed  in  Table  II  which  pertain  to  the  explosive  charge  and 
charge/plate  positioning  must  first  be  specified. 

•  Next,  the  three  soil  parameters  appearing  in  Table  II  (i.e.  seismic  velocity,  soil  fracture 
strain  and  soil  irreversible-compaction  bulk  modulus)  must  be  defined.  However,  as 
implied  by  the  results  in  Figures  3-5,  the  values  of  these  parameters  depend  upon 
both  the  soil  porosity  (a)  and  degree  of  saturation  {(3).  In  other  words,  the  values  of 
these  parameters  can  be  established  only  after  a  and  / 3  are  determined.  The  three  soil 
parameters  mentioned  above  can  then  be  found  using  the  correlation  curves  shown  in 
Figures  3-5.  Thus,  in  order  to  utilize  Figures  3-5  and  calculate  the  expected  impulse 
from  the  resulting  soil  parameter  values,  a  and  f3  must  be  quantified. 

•  To  determine  a  and  (3,  sand/soil  mass  density  (p)  and  the  volumetric  moisture 
content  (w)  of  the  sand/soil  must  be  first  determined  experimentally.  Then,  the 
density,  p  =  (1  -  a)pfuu.compaction  +  a(3pH20,  and  the  volumetric  moisture 
content,  w  =  af3  ,  relations  are  used  to  compute  a  and  j3.  For  convenience, 
functional  relationships  between  a  and  (3  on  one  side  and  p  and  w  on  the  other 
are  determined  fqr  a  standard  value  of  p/uiuompaction  =  2,670  kg/m'J  and 
PH20  =  1,000  kg/m'J  and  depicted,  respectively,  in  Figures  7  and  8;  and 

•  Once  a  and  (3  (and,  in  turn,  c,  Kic,  and  sfr)  are  evaluated,  the  total  impulse  can  be 
computed  from  equation  (2)  as: 
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Figure  7. 

A  correlation  between  soil 
porosity  and  its  density 
and  moisture  content 


Soil  Density,  kg/m3 


Figure  8. 

A  correlation  between  soil 
saturation-level  and  its 
density  and  moisture 
content 


Soil  Density,  kg/m3 
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5.  Summary  and  conclusions 

Based  on  the  results  obtained  in  the  present  work,  the  following  main  summary  points 

and  conclusions  can  be  made: 

•  A  detailed  overview  and  critical  assessment  is  provided  of  the  dimensional  analysis 
proposed  by  Wenzel  and  Hennessey  (1972)  for  sizing  military  vehicle  hulls  required 
to  survive  different  levels  of  blast  loads  associated  with  detonation  of  SBEs. 

•  Using  basic  physical-level  understanding  of  the  key  phenomena  associated  with 
SBE  detonation  and  the  interaction  of  detonation  products  with  the  surrounding 
soil  as  well  as  the  interactions  of  detonation  products  and  soil  ejecta  with  the 
target  plate  (a  surrogate  for  the  vehicle  hull),  it  was  found  that  the  original 
analysis  of  Wenzel  and  Hennessey  (1972)  can  be  greatly  simplified  if  it  is  based 
on  the  physical  phenomena/processes  which  have  a  first-order  effect. 

•  A  new  set  of  non-dimensional  FI  terms  is  introduced  which  not  only  accounts  for 
the  effects  of  target-plate/explosive  charge  geometries,  materials  used, 
deployment  strategies,  and  also  include  critical  effects  associated  with  the  soil 
type  and  structural-condition/consistency. 

•  A  new  functional  relationship  is  introduced  which  relates  the  blast  impulse  to  the 
geometrical,  material,  and  deployment  aspects  of  the  explosive  charge  and  target 
plate  as  well  as  of  the  soil  into  which  the  explosive  is  buried. 

•  A  set  of  blast-impulse  experimental  data  is  finally  used  to  parameterize  this 
relation  and  judge  its  appropriateness. 
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